Q value and half-life of double-electron capture in 184 Os 
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184 Os has been excluded as a promising candidate for the search of neutrinoless double-electron 
capture. High-precision mass measurements with the Penning-trap mass spectrometer TRIGA- 
TRAP result in a marginal resonant enhancement with A = -8.89(58) keV excess energy to the 
1322.152(22) keV 0+ excited state in 184 W. State-of-the-art energy density functional calculations 
are applied for the evaluation of the nuclear matrix elements to the excited states predicting a 
strong suppression due to the large deformation of mother and daughter states. The half-life of the 
transition exceeds T 1/ / 2 ( 184 0s) > 1.3 x 10 29 y for an effective neutrino mass of 1 eV. 



PACS numbers: 07.75. +h, 14.60. Lm, 21.60.Jz, 23.40.-s 

The observation of neutrinoless double-beta transitions 
would reveal physics beyond the Standard Model, as 
it would establish neutrinos to be Majorana particles, 
which implies a violation of the lepton number conserva- 
tion. Experiments searching for these transitions have fo- 
cused on the detection of neutrinoless double-beta decay 
(Of/3/3) rather than neutrinoless double-electron capture 
(Ovee). One reason among others is in general the sig- 
nificantly shorter half-life of the CV/3/3 process. However, 
in the case of neutrinoless double-electron capture, the 
transition is expected to be resonantly enhanced if the 
initial and the final state of the transition are degenerate 
in energy [1—3]. 

In this work, we investigate neutrinoless double- 
electron capture in 184 0s to various nuclear states. In 
particular, the transition to the nuclear excited + state 
with an energy of 1322.152(22) keV by a capture of 
two K-electrons was tagged as partially resonantly en- 
hanced [1]. The capture of two K-electrons, the low spin 
of the nuclear excited state, and the large nuclear ra- 
dius could result in a reasonably short half-life compared 
to other transitions even without being fully resonantly 
enhanced [1]. 

The capture rate of a Ovee transition is given by [1, 3] 
(XG F \V ud \) 4 , ,2| M ££|2p ( r 2^ \ 

(i) 

where A = gA/gv = 1.2701(25) [4] is the ratio of the 
axial vector and the vector coupling constant, Gf the 
Fermi coupling constant, \V u d\ the ud element of the 



Cabibbo-Kobayashi-Maskawa (CKM) matrix, R the nu- 
clear radius, \m ee \ the effective electron neutrino mass, 
|M"| the nuclear matrix element (NME) of the tran- 
sition that consists dominantly of Fermi and Gamow- 
Teller components [5], P ££ = I^mJ l^f/al the absolute 
square of the two electron wave functions at the nu- 
cleus, T 2 h the width of the two-electron hole state in 
the daughter nucleus, and A the excess energy. In or- 
der to estimate the capture rate, we have performed cal- 
culations of the NME by state-of-the-art energy density 
functional (EDF) calculations. These methods have been 
successfully applied to nuclear structure calculations (see 
ref. [6] for a review) and they have been recently imple- 
mented in calculations of neutrinoless double-beta decay 
and double-electron capture NMEs [5, 7]. In order to 
derive A = Q ee — — B211, we have experimentally de- 
termined the double-electron capture Q value Qee > which 
is the difference of the atomic masses of the mother and 
daughter nuclides. E~ t denotes the nuclear excitation en- 
ergy of the final state, and B 2 h is the excitation energy of 
the two-electron hole state in the daughter nucleus. Typ- 
ically, the uncertainty of A is limited by the uncertainty 
of Q ee . In the case of double-electron capture in 184 0s, 
the Q value of Q ee = 1451.2(1.6) keV [8] was only ob- 
tained indirectly from the masses of the other Os and W 
isotopes via (n,7)-reactions. Q values taken from mass 
values in Ref. [8] can be rather inaccurate as the mea- 
surements in the references [9, 10] have revealed. Thus, 
for all potentially resonantly enhanced transitions, their 
Q values must be measured directly with Penning traps, 
which are presently the tool of choice [11]. 
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The resonance condition of several promising transi- 
tions have recently been investigated with Penning traps 
[9, 10, 12-16], and a strong resonant enhancement was 
found for two nuclides, 152 Gd [12] and 156 Dy [13]. 

The half-life for double-electron capture in 184 0s has 
recently been estimated in ref. [1] using an empirical 
value for the NME, which corresponds to the maximum 
value obtained for medium-mass nuclei assuming spheri- 
cal shapes for initial and final states [17]. However, mid- 
shell nuclei like 184 Os and 184 W are known to be signif- 
icantly deformed [18]. As NMEs have been shown to be 
very sensitive to deformation [5, 7], we have calculated 
the NME using a method that consistently describes the 
deformation of the nuclear states involved. We have fur- 
ther recalculated the electron wave functions and atomic 
excitation energies in the daughter nuclide. Hence, we 
report about improvements of three crucial input param- 
eters to the resonance enhancement factor, which allows 
us to check if the 184 Os and 184 W pair is a viable candi- 
date to observe neutrinoless double-electron capture. 

The Q-value measurement reported here was per- 
formed with the double-Penning-trap mass spectrome- 
ter TPJGA-TRAP [19]. The Q value is determined 
by measuring the ratio r of the cyclotron frequencies 
v c = qB/(2wM) of 184 0s+ to 184 W+ ions stored in a 
Penning trap with a magnetic held B: 

Q ee /c 2 = (M( 484 0s) - m e ) (l - ^rS^j) , (2) 

where m e denotes the electron mass. This method pro- 
vides a direct measurement of the Q value. Note that the 
absolute atomic masses M of these two nuclides cannot 
be derived from such a measurement. However, the ab- 
solute masses of the two nuclides were determined within 
this work by measuring the ratio r' = h> c {Ci§) / v c of the 
cyclotron frequencies of the carbon cluster ion 12 C^ 5 and 
the investigated nuclide, in order to obtain the mass val- 
ues via M = (Mc is — m e )r' + m e . 

A laser ablation ion source as described in ref. [20] was 
used to produce 184 0s and 184 W ions from solid sam- 
ples of osmium and tungsten. Due to the low natural 
abundance of 184 0s (0.02%), an enriched sample with 
1.5% abundance of 184 0s was used for the ion produc- 
tion. The ion pulses from the ion source were captured 
in the cylindrical Penning trap, called purification trap. 
There, the captured ion pulse was cleaned from unwanted 
ion species using a mass-selective buffer-gas cooling tech- 
nique [21] with mass resolving power of about 3xl0 4 . 
By using two different targets for 184 W and 184 0s, the 
presence of isobaric contaminants in the purification trap 
was avoided. A cooled monoisotopic ion bunch was trans- 
ferred into the hyperbolical Penning trap, called precision 
trap, where the actual mass measurement took place. 
The cyclotron frequency of the stored ions was measured 
via the time-of-flight ion-cyclotron-resonance (TOF-ICR) 
method with a Ramsey excitation scheme [22] using two 



excitation pulses of 200 ms and a waiting time of 1600 
ms in between. A typical cyclotron resonance for 184 Os + 
is shown in Fig. 1. The cyclotron frequency has been ob- 
tained from a fit of the theoretical line shape [23] to the 
experimental data. 
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FIG. 1. Time-of-flight ion cyclotron resonance of 1 Os + with 
about 450 ions. The data points show the mean time-of-flight 
values as function of the excitation frequency. The solid line 
is a fit of the theoretical line shape [23] to the data points. 

In order to determine the Q value according to Eq. (2), 
the cyclotron-frequency ratio r of 184 0s + to 184 W + was 
determined in alternating measurements of the individual 
cyclotron frequencies (see Fig. 2). In the evaluation pro- 
cedure for cyclotron frequency ratios, systematic uncer- 
tainties as described in detail in ref. [24] were considered. 
The systematic shift of the frequency ratio depending on 
the mass difference Am between the two ion species used 
in the measurement is negligible (e m /r < 10~ u ) in case 
of nuclides with the same mass number A, as in a Q- value 
measurement of a Ovee process. For the measurement of 
the individual masses, this mass-dependent effect with a 
mag nitude of e m = 2.2(0.2) x 10" 9 Am/u was considered. 
It arises from a possible slight misalignment between the 
electric and magnetic fields of the Penning trap [24]. 

The result of the cyclotron-frequency ratio measure- 
ment of 184 Os + to 184 W + and the Q value derived with 



TABLE I. Results of the mass measurements. For each ion 
of interest the reference ion species is 12 C 1 " 5 . The cyclotron- 
frequency ratio r', the mass excess ME = M - A u, where u 
is the unified atomic mass unit, and the mass excess MEut 
listed in Ref. [8] are given. 

Ion r[ ME (keV) ME nt (keV) 

184 Os+ 1.0219583764(67) -44251.47(1.13) -44256.60 (1.34) 
184 W+ 1.0219496960(55) -45705.40(0.94) -45707.54 (0.87) 
12 Cje 1.0666668659(70) -0.62(1.25) (0) 
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FIG. 2. The data points show the 32 cyclotron-frequency ratio 
measurements recorded to determine the double-electron cap- 
ture Q value of 184 0s. The thick line represents the weighted- 
mean value and the gray area its 1 a-uncertainty. 

Eq. (2) are 0.9999915163(38) and 1453.68(0.58) keV, re- 
spectively. The uncertainty of the literature Q value was 
improved by a factor of 3. A slight deviation below 1.5 a 
from our result to the value in Ref. [8] (Q ee — 1451.2(1.6) 
keV) was found. 

The atomic mass excess values obtained from the 
cyclotron-frequency ratio measurements in this work are 
listed in Tab. I. Our measurement is a direct compar- 
ison to the atomic mass standard, unlike the values 
from Ref. [8], which show a deviation of 5.13(1.75) keV 
and 2.14(1.28) keV for 184 0s and 184 W, respectively. 
The Q value extracted from our values ME( 184 0s) and 
ME( 184 W) listed in Tab. I is 1453.93(1.47) keV and 
agrees well with our direct Q- value measurement. 

In order to determine the excess energy A to the ex- 
cited 0+ state with 1322.152(22) keV energy in 184 W 
[25], the excitation energy of the two-electron hole state 
£?2h and its width I^h have been calculated. The 
value of B 2 h was obtained by applying the Dirac-Fock 
method [26], where frequency-dependent Breit interac- 
tion, quantum electrodynamics (QED) and correlation 
corrections were included. The calculations were per- 
formed for the Fermi model of the nuclear charge distri- 
bution (with r RMS =5.3670 fm for 184 W [27]), resulting 
in B 2h = 140.418(12) keV. The value for the width of the 
autoionizing state of 184 W with two K-shell holes is 
approximated by using the double value of the width of 
the single K-level hole [1]. The single K-level hole width 
is taken from Ref. [28]. We estimate the uncertainty of 
T2h by using the most conservative value in Ref. [28] and 
increase it by 2 eV to account for the uncertainty due to 
the hole-correlation effect suppressed approximately by 
a factor 1/Z. The contribution of the width of the ex- 
cited nuclear state has been neglected. This approach 



resulted in = 80(10) eV. Prior to this work, A was 

— 11.3(1.6) keV using Ref. [8] for the two mass values. 
The uncertainty given is obtained by adding the individ- 
ual uncertainties of 1.3 keV and 0.9 keV of 184 0s and 
184 W, respectively. Using our measured value, A be- 
comes -8.89(58) keV. 

The evaluation of the NME with the EDF approach 
was performed with the Gogny DIM functional [29]. In 
our EDF approach the various initial and final nuclear 
states are given as superpositions of particle number and 
angular momentum projected Hartree-Fock-Bogoliubov 
states |\I/^(/3)) with different quadrupole deformation (3 
(P < 0, /3 = and (3 > represent oblate, spherical and 
prolate shapes respectively) [7] : 

i J J/ r ) = Efv/C 8 )i*i//C8))' ® 

where are the coefficients of the linear combi- 

nation which are obtained by solving the corresponding 
Hill- Wheeler-Griffin (HWG) equations and a labels the 
different states for a given angular momentum / [31]. 
From these coefficients, the so-called collective wave func- 
tions - the amplitude of the states for a given deformation 

- can be extracted [31]. Fig. 3 shows these amplitudes in 
the three lowest + states in 184 W and the ground state 
in 184 0s computed using the Gogny DIM functional. All 
these states are well prolate deformed. The two ground 
states have a similar distribution around /? ~ 0.25; the 
first excited + state of 184 W correspond to a vibration of 
the deformed ground state while the 0jj~ state has a larger 
prolate deformation j3' ~ 0.35. These amplitudes will af- 
fect significantly the final value of the NMEs as we will 
analyze below. In Table II we show the excitation energy 
of the first excited 2 + state E(2f) and its quadrupole 
transition probability B(E2) to the ground state. We 
find larger excitation energies of the first 2f state and 
reduced quadrupole transition probabilities compared to 
the experimental values. Tab. II lists also the energies 
of the two lowest excited + states in 184 W. While the 
excitation energy of the 0^ state has a fair agreement 
with data, the 0jj~ excitation energy is too high, pointing 
to the need of including additional degrees of freedom. 
For example, the inclusion of pairing vibrations coupled 
to the quadrupole degree of freedom would lower the en- 
ergies of the 0^,0 J states, bringing them closer to the 
experimental values. However, the quadrupole deforma- 
tion of these states would not be very much affected [30] . 

The NMEs from the 184 0s ground state to the three 
different + states in 184 W are obtained in this frame- 
work by folding the deformation dependent wave func- 
tions with the Fermi and Gamow- Teller matrix elements 
expressed in the deformed basis [7]: 

M«(0+) = (O+f^Mjf/^loXt), (4) 
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E(0+) E(0+) 


£(2+)[ 184 0s] £(2+)[ 184 W] B(£2)[ 184 0s] B(S2)[ 184 W] M"(0+) 


M«(0+) M«(0+) 


M"(Sph.) 


Expt. [32] 
EDF DIM 


1.002 1.322 
1.135 3.169 


0.120 0.111 99.6 119.8 

0.172 0.155 135.0 146.6 0.631 


0.504 0.163 


8.177 



TABLE II. Excitation energies (in MeV) of the first two + states for 184 W and of the 2+ state and its quadrupole transition 
probability (in Weisskopf units) for 184 Os and 184 W are compared with EDF calculations based on the Gogny DIM functional. 
The values of the NMEs, M", to the three + states in 184 W and the NME obtained assuming spherical shapes (/3 = 0) are 
also included. The experimental resonant state is written in boldface. 
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FIG. 3. (Color online) Deformation amplitudes for the three 
first + excited states of 184 W (a) and the ground state of 
184 0s (b) calculated with EDF Gogny DIM. (c) Intensity of 
the Fermi component of the Ovee NME as a function of the 
quadrupole deformations (3 and /3' of the 184 Os and 184 W 
nuclei, respectively. 



where MJy GT are two-body Fermi and Gamow- Teller 
transition operators. The results are listed in Table II. 
We obtain that the Gamow- Teller component contributes 
about 75% to the NME. For comparison we have also 
listed the NME obtained if we hypothetically assume 
spherical shapes for the initial and final nuclear states 
(i.e. /? = f5' = 0), obtaining values close to those of 
ref . [1] . However, if the deformation of the nuclear states 
is properly accounted for, the NMEs are reduced by 
about an order of magnitude for the ground state tran- 
sitions and even more for the transitions to the excited 
states. To understand this behavior we have plotted in 
Fig. 3(c) the Fermi component of the NME as a func- 
tion of the quadrupole deformations of the 184 0s (/?) and 
184 W nuclei (the Gamow- Teller component shows 
the same general behavior): 



:^Bi o w (/3oi*fe o w (^o)(^ o 0s (/3)i*fe o o s (/ 3 )) 



(5) 



M 



F/GT 



We observe two general trends: a) for identical deforma- 
tions of initial and final states (i.e. /? = /?') the NME 
decreases with increasing quadrupole deformation; b) for 
fixed deformation j3 of the initial state the NME decreases 
with increasing difference between the deformations of 
initial and final states. We find significant and similar 
quadrupole deformations for both ground states, consis- 
tent with the large B(E2) values observed experimen- 
tally. This explains why our NME values are so much 
reduced compared to the assumed value in ref. [1]. Al- 
though the first excited + state in 184 W has a similar de- 
formation as the ground state, the corresponding NME is 
smaller. This is caused by the fact that the wave function 
of the first excited state has to be orthogonal to the one of 
the ground state implying a nodal structure which leads 
to cancellations when calculating the NME. In fact, we 
could interpret this state as a vibrational excitation of the 
ground state. The second excited + state is predicted by 
our calculation to have a noticeably larger deformation 
than the ground state. This results in a smaller NME 
than for similar deformations as obtained for the ground 
state due to the trends observed in Fig. 3(c). Even if our 
calculation fails to predict the correct deformation of the 
0$ state, the NME to this state should be smaller than 
the one to the ground state. If its deformation is differ- 
ent than the one of the ground state, the general trend 
b) predicts a smaller NME. If the deformation is similar, 
there should be strong cancellations in the NME caused 
by the orthogonality to the two lower + states. In this 
way, the NME for the capture to the ground state can be 
interpreted as an upper limit. 

The half-life of the resonant neutrinoless double- 
electron capture in 184 0s is related to Eq. (1) by Ty2 — 
In2/A ee . Using R = 5.382 fin from [27], and P ee = 
1.196 1 7 x 10 33 eV 6 (in natural units) determined by the 
one-configuration Dirac-Fock method for an extended nu- 
cleus [26] and the average of the two values for M e€ (0^") 
in Tab. II, we obtain TWa = 2.0 x 10 30 y for an effective 
neutrino mass of 1 eV. If we take the ground state NME 
value as an upper limit (see above) , we find a lower limit 
for the half-life of T 1/2 = 1.3 x 10 29 y. 

In conclusion, we have investigated the resonance con- 
dition and the NME of double-electron capture in 184 0s 
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to three 0+ states in 184 W. The Q value was determined 
in a direct measurement of the mass difference of 184 Os 
to 184 W to be 1453.68(0.58) keV with a factor of 2.7 
smaller uncertainty than previously reported. We have 
computed the NME in a microscopic approach that con- 
sistently considers deformation degrees of freedom and 
found that initial and final nuclear states are significantly 
deformed. This deformation results in a suppression of 
the NME by a factor of about 10 compared to a calcu- 
lation that is restricted to the spherical configuration. 
Combining our various improvements, a lower limit of 
the half- life for the Ovee process of 1.3 x 10 29 y to the ex- 
cited 0+ state with 1322.152(22) keV excitation energy m 
184 W, normalized to an effective electron neutrino mass 
of 1 eV, has been obtained. Thus, 184 0s can be excluded 
as a promising candidate for the search for neutrinoless 
double-electron capture. 
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